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Abstract 

Behavior of doped fermions in Z2 gauge theories for the quantum dimer and eight-vertex models is 
studied. Fermions carry both charge and spin degrees of freedom. In the confinement phase of the Z2 
gauge theories, these internal symmetries are spontaneously broken and a superconducting or Neel 
state appears as the groundstate of the doped fermions. On the other hand in the deconfinement- 
topologically-ordered state, all symmetries are respected. ^^From the view point of the quantum dimer 
and eight-vertex models, this result indicates interplay of the phase structure of the doped fermions 
and that of the background dimer or the eight-vertex model. At the quantum phase transitions in 
these systems, structure of the doped fermion's groundstate and also that of the background dimer or 
eight-vertex groundstate both change. Translational symmetry breaking induces a superconducting 
or antiferromagnetically ordered state of the doped fermions. 

1 Introduction 

Study of the strongly-correlated electron systems is one of the most interesting problems in the con- 
densed matter physics. Most of the strongly-correlated electron materials like the high-temperature(T) 
cuprates, the heavy fermions, etc. have various phases like the superconducting state, antiferromag- 
netic(AF) Neel state, etc. in very close parameter regions. This fact indicates that in those materials 
the charge and spin degrees of freedom of electrons correlate and interplay with each other and as 
a result the phase transitions occur. Then it is important to understand quantum phase transi- 
tions(QPTs) which are difficult to be described by the ordinary Ginzburg-Landau theory in terms 
of a single order parameter. Theoretical studies on the QPTs are welcomepP but at present there is 
no unified and transparent understanding of the QPTs which connect various phases with different 
order parameters. 
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In this paper we shall investigate the QPTs in which various order parameters interplay with 
each other. In particular we are interested in the quantum dimer and eight-vertex models with the 
doped fermions. Roughly speaking, the "background" dimer state can be regarded as f-electrons 
and the doped fermions are electrons in a conduction band from the view point of the heavy fermion 
materials. As the dimer states can be described by a gauge theory as it is well-known 0, the 
fermion doped system under study is described by a Z2 Ising gauge theory(IGT) coupled with the 
fermions. The doped fermion has hoth the charge and spin degrees of freedom. Then it is interesting 
to see if there exists an interplay between background dimer states like valence-bond crystal or liquid 
and the order of charge and/or spin degrees of freedom of the doped fermions. 

In the elementary particle physics, it is well know that confinement of quarks induces spontaneous 
chiral symmetry breaking, i.e., phase transition of confinement and that of the chiral-symmetry 
breaking occur at the same time. Then one can expect that similar interplay of different phase 
transitions occurs in the present model. We shall study this interesting possibility in this paper. 

The present paper is organized as follows. In Sec. 2, we shall study the IGT coupled with fermions. 
In particular behavior of the charge and spin degrees of freedom of the fermions is investigated 
in strong-coupling regions. We show that in the even IGT the superconducting state is realized, 
whereas in the odd IGT the Neel state appears as the groundstate. On the other hand in the 
weak-coupling region, any symmetry breaking does not occur. In Sec. 3, we apply the above result 
to the quantum dimer model on the square lattice. In the valence-bond crystal phase, a dimer or 
the Neel like state of the fermion's spin degrees of freedom appears. In Sec. 4, we consider the doped 
quantum eight-vertex(q8v) model. In the ordered state of the q8v model, the doped fermions are 
in the superconducting phase, whereas in the disordered phase low-energy excitations are weakly 
interacting fermions. Section 5 is devoted for conclusion. 

2 Symmetry breaking in the IGT with fermions 

In this section, we shall first review the IGT and study behavior of fermions coupled with IG 
variables|3|. We use the Hamiltonian formalism on d-dimensional hypercubic spatial lattice. 

The IG variables are defined on links of the lattice {x, £) where x denotes site and the direction 
index £ — 1, • • • , d. Obviously {x, £) = {x + £, —£) where £ is the unit vector of the £-th direction. 
Quantum state of the IGT is specified by the quantum state at each link which is a linear combination 
of "up" and "down" states; ci| f) +C2I [) = (ci,C2). Then the IG operator is given by cr|(a:), where 
<y'^{x) is the Pauli spin matrix acting on the gauge-field configuration (ci,C2) at {x,£). Similarly 
conjugate variable, i.e., the electric field operator, is given by aj{x). Hamiltonian of the IGT is 
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composed of the electric and magnetic terms, 

-r^ai-At^dVW, (1) 

link pi 

where (link) and {pi) denote Hnks and plaquettes, respectively. F and k in Ti.^ are parameters and 
they are related with the gauge coupling g a,s T (x and k cx l/g^- The above Hamiltonian is 
supplemented by the physical-state condition. In the even(E) IGT, 

= Yl (ji; G.,\phys) = \phys), (2) 

link^x 

where {link e x) denotes the links emanating from the site x. The phase structure of the EIGT is 
well known. For d — 1, there is no phase transition and only confinement phase exists. On the other 
hand for d = 2, 3, there is a phase transition from confinement to Higgs phase as T/k decreases. 
The odd(O) IGT is defined by the same Hamiltonian (Q, whereas the local constraint is changed 

to, 

Gx\phys) ^ -\phys). (3) 

It is easily verified that the above constraints ©and ^ commute with He- The OIGT is closely 
related with the dimer model as we explain later on^. Detailed phase structure of the OIGT is not 
known yet. 

Let us introduce fermions ip°'{x) {a — 1,2, ■ ■ ■ , S) and couple them to the IG field crf (a;). Fermion 
part of the Hamiltonian is given as 

n^ = -t J2 ^i{x)ai{x)r{x+i)+ii.c+MY,^ii^)r{x), (4) 

link.a x 

where t is the hopping parameter and the "chemical potential" M controls the fermion density 
(ip'^ip). The total Hamiltonian of the system is given by 

Ti-T = 'He + H^. (5) 

By the fermion coupling, the physical-state condition changes to 

. e'" So ("^'^^^^ \phys) = \phys) , for EIGT (6) 
■ e"I^o'/''(^)'/'(^)|p/iys) = -\phys), for OIGT. (7) 

Let us consider the EIGT first. Effect of the coupled fermions to the phase structure will be 
discussed in the following section. It is shown there that the coupling of the fermions enhances 
the deconfinement phase but there still exists the confinement-deconfinement phase transition. In 
the deconfinement phase k/F ^ 1, fluctuation of the gauge field (j'j{x) is small and we can simply 
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put <jf{x) ^ 1 (up to irrelevant gauge transformation). The ferniions move almost freely and no 
spontaneous symmetry breaking(SSB) of the internal symmetries occurs. On the other hand in the 
confinement phase F/k 3> 1, one can expect that a SSB occurs as in the most of the gauge theories. 
A good example is the SSB of chiral symmetry in the quantum chromodynamics. 

In order to investigate the SSB in the present gauge model, we shall derive a low-energy effective 
model of fermions in the confinement limit F/k 3> In the leading order of F, the lowest-energy 
states are specified as crKx) = — 1 for all links and the fermionic sector is highly degenerate. In order 
to derive an effective model for the fermionic sector, we introduce the projection operator V to the 
subspace \iIj)gs satisfying <j}{x)\iP)gs = —Was for all hnks, and we put Q= l — V. We start with 
the eigenvalue problem for the original gauge-fermion system, 

nTW = E\ij). (8) 

By using V and Q, we can derive the following equations, 

{VHtV - EV) IV) = -VHtQ W 
{Q-HtQ - EQ) = -Q-HtV IV) 

Q\'<P) = {E-QnTQ)-^QnTV\^) (9) 
^From Eqs.®, we have an eigenvalue equation in the subspace T'lV') — I'0)g5'j 

VUtV + PHTQQiE ~ QUtQ^^QUtV] |V) = EV |V) ■ (10) 
Then an effective Hamiltonian He{E), which acts on the subspace \4')gSi is given by, 

He{E) = VHTr + VHTQQ{E-QHTQr^QHTP- (11) 
Explicit form of He{E) is obtained in the leading order of 1/F as follows, 

He{E) ^ -^5][vi(x)V"(^ + ^>^(x)V''(^ + ^')" 

x,t 
x,i 

-^Y.[^'^o.{x+l)r{x)4{x)i^'{x + i)\ (12) 

x,l 

where we have used 

T^W^{x)]Q{^l{x')}V = V^-,-'^- (13) 



4 



He (E) acts on the fermion sector of |^/')gs and resolves its degeneracy. 

Let us comment on the above low-energy effective model for the strongly coupled gauge theory. 
^From the Hubbard model at the half filling, an antiferromagnetic(AF) Heisenberg model is derived 
as a result of the strong on-site repulsion. The above effective Hamiltonian Hc,{E) H12I) also contains 
the AF magnetic interactions as we see shortly. Origin of the AF magnetic interactions is not the 
on-site repulsion but the strong attractive gauge force. One may expect that this attractive force 
may induce other SSBs besides the long-range AF order, like the superconductivity. This is the case 
as we see in later discussion. 

The fermion fields satisfy the following canonical anticommutation relations 

{r{x),i^ny)}^Q. (14) 
Then it is useful to define the "spin" and "charge operators", Q(a;)'s and /'(a;)'s, 



(15) 



W = \ rix)] - i<5|Q(x) = 4{x)r{x) - \5f, - \5fMx) (16) 

Pcp{x) = \[4'i{x),4{x)\=^ki{x)^l^l{x) (17) 

Po.p{x) = \[r{x),^^{x)\=r{x)i'^{x). (18) 

Commutation relations of Q(x)'s and P(a;)'s are readily obtained from Eas. H14|l . 

[Q{x),Q{x)] = [Pa/3 (a;), Pa/3 (x)] = \P o^pix) ,P c.i3{x)] = (19) 

[Q1{x),Qj{x)] = S2Q}{x) - SjQfix) (20) 

[Q(x), Pp-,{x)] = 2Pp^ = -2P^0 (21) 

[Q{x),Pp^{x)]^-~2Pp-,^2P.,p (22) 

[PMx),Pys{x)] = S;Qiix) - 8lQ%(x) - 5iQ}{x) - S^QUx). (23) 

In terms of Q(x)'s and P(a;)'s, He{E) can be written as follows, 

3+2 +2 _ _ 

He{E) = - — V+—Y,[PMx)Pc.0{x + e)+Pc.f3ix)P^0ix + i)] 



^ E[2Q^(x)g^(.T + i) + Qix)Qix + i)], (24) 

x,i 



where V is the total number of sites of the lattice. 
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Let us consider the simplest case in which the parameter a takes two values, i.e., a =t, i-^ Then 
the states at each site are explicitly given by 

hi), V'll-l), ^jV-Il-l) (25) 

where | — 1) is the empty state, — 1) = and Qj — 1) = — | — 1). 

In the EIGT, only two out of four states | — 1) and i/'j-i/'l | — 1) are the physical state for aj (x) = — 1 
and by Eq.®. Then in the Hamiltonian He{E), the terms Q'^{x)Q^{x + £) in do not contribute. 
The groundstate of He{E) can be written as 

1 5) = n ["(2^)1 - 1).^ + vix),p\{x)i;lix)\ - 1},] , (26) 

X 

where u{x) and v{x) are complex constants which satisfy \u(x)\'^ + \v{x)\'^ = 1. Obviously, the above 
state \S) is nothing but the superconducting state. It is straightforward to calculate the expectation 
value of the energy for the state H26(l . and the lowest-energy state is obtained as u{x) — u and 
v{x) = e{x)v where t{x) = 1(— 1) for the even site(the odd site). Then the energy Es is obtained as 

Es - + ^ 5]h4..V + i-u' + vYl (27) 

On the other hand the density of fermions is 

'J2^ir)^{{Q + i)) = '2v'- (28) 



Therefore at the half filling, = -y^ = 1/2 and Es = -^^V. 

It is not so difficult to show that the same effective Hamiltonian H24() is derived in the OIGT at 
half- filled case. In this case, exactly one fermion resides at each site and the state (jI{x) = — 1 is 
realized for all links. In the single-particle sector of the fermions at each site {-(Ajl — I),''/'}! ~ 1}}, 
only the last two terms in He{E) (|24|l operate and the model reduces to the AF Heisenberg model. 
At present it is established that the Neel state with the long-range AF order is the groundstate for 
d > 2. Therefore the groundstate of the doped fermions in the half-filled OIGT is the AF magnet. 

In the following section, we shall study the phase structure of the doped IGT by the mean-field 
approximation(MFA) before going into detailed investigation on the doped dimer and q8v models. 



3 Phase structure of the doped IGT 

In this section we shall study the effect of the doped fermions to the phase structure by the MFA. 
We first consider the pure EIGT. It is not so difficult to derive MF Hamiltonian Hm from We 

^Please do not confuse the fermion index a =|, J, with that of the Z2 gauge state. 
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formally decompose (Tj{x) into the MF Uq and the "fluctuation" Saf{x) from it, crf{x) = Ua + 5(7^{x). 
The Hamiltonian Tig. of the IGT is written as 



link 



link 



link 



link 



(29) 



where Np and Nl are the number of plaquettes and links of the lattice and for the d-dimcnsional 
hypercubic lattice Nl — -^jziy-ZVp. ^From Ea. (|^ and by neglecting the terms o{^5a\(x)Y^, we 
obtain the MF Hamiltonian 

Uii^vY, a\ (x) + K [SU^NP ~ 2{d - l)C/o' ^ (^)] ' (^0) 

link link 

^From Ea. (|3U|l . it is straightforward to calculate the energy of the groundstate as a function of the 
MF Uq. Here we introduce the gauge coupling constant g. As we stated before T = and k — l/g^. 
The effective potential, i.e, the energy of the groundstate per link, V^piUo) is obtained as 



V^ipiUo) 



3{d~l) 
~^2g^ 



(31) 



In Fig^, we show V^p{Uo) for various values of the gauge coupling g^. Result obviously shows that 
there is a first-order phase transition and at weak coupling ((t|(x)) ^ 0, i.e., the phase corresponding 
to the Higgs phase appears^. At present for the 3D (d — 2) IGT, it is known that there exists a 
second-order confinement-Higgs phase transition whereas in the AD {d — 3) IGT the transition is of 
first-order 

Effect of the doped fermions is investigated easily by the MFA. We simply put crKx) ^ Uq in 
the Hamiltonian 7t!^ Q to obtain a MF Hamiltonian H 



M' 



nlj = -tUo V^(x)V'(2: + i) + H.c. + mJ^^^^M^)- (32) 

link X 

The MF Hamiltonian can be easily diagonalized by the Fourier transformation of the field operators, 
T^ti = f{dp)4'Hp) \ - 2tUo{J2 ^o^^'^) + ^'^1 (33) 



The energy of the groundstate per site is given by 



COS Pi 



(34) 
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where the Fermi momentum pp is determined by the density of the fermions pp, e.g., pp ~ 
{i-n"^ PfY^'^ for d = 3 (and positive Uq). iFrnva H34|l . it is obvious that V^jp{Uo) oc — t|?7o| and 
therefore the fermion doping enhances the condensation of af, i.e., the deconfinement-Higgs phase. 
This resuh by the MFA, i.e., the hnear-f/g dependence of V^^(C/o)i probably over estimates the 
fermion effect. But we expect that the above result of the enhancement of the deconfinement phase 
is qualitatively correct as it is observed in the various gauge systems.^ 
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Figure 1: Effective potential V^jp{U[)) of the pure EIGT with the gauge coupling g = 0.6, ■ • ■ , 1.2. 
The result indicates a first-order phase transition at g ~ 0.9. 



4 Doped quantum dimer model 

In this section, we shall consider quantum dimer model(QDM)[5j which was introduced as a model of 
resonant-valence-bond (RVB) statelH). We consider the two-dimensional square lattice for simplicity 
in the rest of this section. Quantum states of the model consist of closed-packed and hard-core 
dimers. Hamiltonian of the quantum dimer model is given as 

Hdimer = ^ [ - Q^P + pVp , (35) 
plaquette 

where J^p is a flip term defined on each plaquette (see Fig|21) and Vp is the potential term which 
takes unity for flipablc dimer state in the plaquette and otherwise vanishing, q and p are parameters 
■^Improvement of the simple MFA was discussed in various places. See for example, Ref . 171 . 
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of the model. At the Rokhsar and Kivelson(RK) point q ^ p, the groundstate is given by the 
summation of all dimer configurations with equal weight; the valence-bond liquid state. 

I 

Figure 2: Flip of a pair of dimers in a plaquette. Solid lines denote dimer and dots are sites of the 
spatial lattice. 

In the previous section we introduced the OIGT. The OIGT can be regarded as a kind of the 
QDM. In the OIGT description of the QDM, the state (jI{x) = —1 on the link {x,£) corresponds 
to the empty state, whereas that of crj{x) = 1 is regarded as the state occupied by a dimer. Then 
the physical-state condition of the undoped OIGT with the local constraint iPJ requires that one or 
three dimers emanate from each site. For large positive F in the Hamiltonian the space of states 
consists of closed-packed and hard-core dimers on the square lattice. On the other hand, the second 
term of Eq.lQJ works as the flip term of two dimers that reside on a single plaquette (see Fig|2). 

At present it is known that in most of the parameter region the columnar or staggered ordered 
state, which breaks the translational symmetry, is the groundstate of the undoped QDM Hdimer 
on the square lattice^Ulj whereas at the RK point the translationally symmetric state is realized|21 
111! I12| . However by the doping, a liquid state of VB may appear as the groundstate in a finite 
parameter region[T^. 

In the previous section we considered the half- filled case of the doped OIGT. Here we focus 
on the lightly doped QDM assuming that a valence bond crystal like the columnar state is still 
the groundstate. Let us consider a two-fermion doped QDM first. A pair of fermions doped in 
the columnar state reside on the nearest- neighbor (NN) sites and a dimer on the link connecting 
the two sites disappears. This state is the lowest-energy state. Similarly when 2n fermions are 
doped into the QDM, there are various configurations in which n dimers disappear by the doping. 
Effective Hamiltonian of the doped fermions can be derived in the framework of the OIGT and the 
Hamiltonian of the AF Heisenberg model is obtained (up to an irrelevant constant term), 

^d^e. = ^ ^ {2Q''^{x)Qi{x + i) + Q{x)Q{x + 1)) , (36) 

where the summation is performed over the links (x, t) connecting sites of doped fermions. 

For an isolated pair of fermions, it is easily to obtain the wavefunction. In particular in the case 
5 = 2, i.e., a =t, |, they form a spin-singlet state, 

(^1 (x)^| {x + i) - (x)^j (x + i)) I - 1), I - 1)^^,-, (37) 
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and the energy of the spin-singlet pair is estimated from (jlffill as Essp = ^t^- Two fermions 
can be put on two sites which separate more than one lattice spacing if the columnar-staggered- 
mixed configuration of dimcrs is considered (see FiglJJ. However those states have higher energy 
compared with the spin-singlet pair state 1)37(1 regardless of their spin configurations as seen from 
jjdimer ^ In this way, we can expect that lightly fermion-doped dimer system has the groundstate 
with inhomogeneous fermion density, i.e., high and low fermion density regions appear. In the 
high density region, the Neel-like state for the fermion's internal degrees of freedom appears or a 
deconfinement-VB-liquid phase is realized. On the other hand in the low density region, a closed- 
packed and hard-core dimer state appears. Numerical studies are required in order to obtain detailed 
quantitative phase diagram of the doped dimer systems. 



(a) Columnar 



(b) Staggered 























(c) Mixed 1 























(d) Mixed 2 



Figure 3: Configurations of fermions in dimer backgrounds: Heavy lines denote dimer. Up(down) 
arrow denotes spin-up(down) state of the doped fermion. 



5 Doped quantum eight-vertex model as IGT 
5.1 Eight- vertex model: Short review 

In this section we shall investigate the behavior of doped fermions in the quantum eight-vertex(q8v) 
model 14 . Classical 8v model was studied by Baxter very intensively JS]. Dynamical variables are 
"arrows" sitting on each link of the square lattice, and then at each site four arrows meet; i.e.. 
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vertex. At each vertex, the numbers of incoming and outgoing arrows are both restricted to even, 
and therefore there are eight types of vertices, see Fig0| For each vertex, the Bohzmann weight 
Wi (i = 1, 2, • ■ • , 8) is asigned and total Bohzmann weight for a configuration is simply given by the 
product of Wi of each vertex composing that configuration. We choose the weight Wi as follows, 

wi ^ W2 ~ a, UI3 = UI4 — 6, 

^ We = c, wy = wg, ~ d. (38) 

Hereafter, a, - ■ ■ ,d also denote type of vertices. Then the partition function of classical 8v model is 
given as 

Z^sv{a,b,c,d) ^ a''"6"'c''<^d''^ (39) 

conf's 

where Va is the number of the a-type vertex in the configuration, etc. 

Basis of the Hilbert space of the q8v model consists of the configuration space of the classical 
8v model. Different configurations of vertices are orthogonal with each other and the norm of each 
configuration of vertices is normalized as unity. Hamiltonian is composed of a fiip term and a 
potential term. The q8v model can be described by an IGT, which is the subject of the following 
subsection. 



aabbccdd 

Figure 4: Vertex configurations: There are eight vertices in eight-vertex model. Their Boltzmann 
weight is given by a, ■ • ■ , d. 

5.2 IGT for q8v model 

As explained in the previous subsection, state of the q8v model is specified by assigning arrow for 
each link of the two-dimensional square lattice. We call the states of the up and right arrow spin-up 
state whereas the states of the down and left arrow spin-down state. In order to construct physical 
operators, let us introduce the Pauli spin matrices for each link {x,£) {£ — 1,2). Let the above 
up-spin and down-spin states on the link (x^i) be eigenstates of the operator cr|(a;) with eigenvalue 
±1, respectively. ^From the restriction on the configurations of the classical 8v model, the physical 
state of the q8v model must be satisfied the following local constraint, 

al{x)(Tl{x — l)al{x)al{x — 2)\phys) = \phys), for all x. (40) 
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The constraint (|40|l is nothing but the gauge-invariant condition in the IGT. Flip operation works 
on arrows on four hnks forming a plaquette, and then the flip term of the Hamiltonian is given as 

HfUp = -5](T3(a:)a3(.T + 2)alix)alix + l). (41) 
pi 

It is obvious that the above Hamiltonian H41|) is the magnetic term of the IGT and commutes with 
the constraint H40(l . 

If there are no terms in the Hamiltonian besides the flip term H41|) . the groundstate is the 
summation of all the states in the configuration space with equal weight, a ^ b = c = d. This state 
is nothing but Kitaev's state which may play an important role in quantimr tori codesTG, . In order 
to give nontrivial weight for states in the q8v model, we shall introduce a potential term. To this 
end we consider the following operators that distinguish the verices a, • ■ • , d, 

Saix) = \{cTlix) + alix) + al{x - i) + alix - 2)) , 

Sb{x) = ^[al{x)~<jl{x)+al{x-i)-al{x-2)), 

Sc{x) = l{'yl{x) ~ alix) - alix - 1) + al{x - 2)) , 

Sd{x) ^1{<tI{x) + alix) - alix -I)- a^x - 2)). (42) 

Operation of the above four operators is summarized in Tabled ^From the Tablet, it is obvious 
that all eight vertices are distinguished by Sa^ • • ■ , Sd- Projection operators can be easily obtained, 

V.^Sf, i^a,b,c,d, (43) 

and after some calculation 

Va+Vb + Vc + Vd^ I. (44) 



Table 1: Relation between operator Si and eight- vertex configurations 
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By using the projection operators Vi (i — a,b,c,d), we can construct a potential term which 
assigns weight for each pure state just like the partition function of the classical 8v model To 
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this end, let us see how vertices interchange with each other under the plaquette flip. For example 
the d-type vertex on left-down or right-up site(even site) of a flip plaquette changes to the a-type 
vertex. On the other hand, the d-type vertex on right-down or left-up site (odd site) changes to the 
6- type vertex. These are summarized in Table|21 We can give the potential term|14|. 

^-^ \a c a / 

X 

X {-Va{x + i) + '{vb{x + i) + -Vc{x + i) + ^Vdix + i)) 

\a c d / 

X (-Va{x + 2) + ^Vb{x + 2) + -P,{x + 2) + ^Vd{x + 2)) 
\a c d / 

X (-Va{x + i + 2) + ^Vbix + i + 2) + -Vcix + i + 2) + ^Pdix + 1 + 2)). (45) 
\a b c a / 

The total Hamiltonian of the q8v model is then given by 

Hqav = Hfup + Hv (46) 

It is not so difficult to show that weight of the configuration with Va a-type vertices etc. is given by 
a^'^Jf^c^'^d"'^ in the groundstate wavefunction \GS{a^ b, c,d)). Then the partition function of the q8v 
model at T = is related with the c8v model as follows, 

Zqsv{a, b, c, d) = Zcs,v{a^, b^,c^,d^). (47) 

Phase structure of the classical 8v model is well-known, and from the relation 1)471) it gives that of the 
q8v model at T = 0. FigElshows the phase diagram of the qSv modeljElE]- We put a = 6 = 1 for 
simplicity. The region / (//), which is given by > d^ + 2 (d^ > + 2), is an ordered phase where 
the c-type vertices (d-type vertices) dominate. On the other hand the region ///, — 2 < — d^ < 2 
is a disordered phase and various states including all vertices appear in the wavefunction of the 
groundstate^. In particular at Kitaev's point a = 6 = c=:d==l,all configurations appear with the 
equal weight. 

Let us consider fermion doping in the qSv model. In the region ///(the disordered phase), doped 
fcrmions move almost freely and no SSB occurs. On the other hand in the ordered phases / and //, 
strong interactions work on doped fermions and we expect that some kind of SSB of the internal 
symmetry of fermions takes place there. This is the subject of the following subsection. 

^In certain limits of the parameter space like = d? ^ oo, only c and/or d vertices appear. See later discussion. 
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Table 2: Relation between flip operation and eight-vertex configurations: Cross symbol and circled- 
cross symbol denote a-type and 6-type vertices, respectively. Similarly, dot symbol is c-type, and 
triangle symbol is d-type vertices. Relation between types of vertices and symbols is summarized in 
the left table. The right table shows how vertices interchange with each other under plaquette-flip 
operation. 
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Figure 5: Phase Diagram of q8v model. Regions / and // are ordered-solid phase whereas region 
/// corresponds to disordered-hquid phase. Doped fermions are in supercondcting phase in region 
/ and //. Region /// are normal phase of doped fermions. 

5.3 The limit —>■ oo with finite cP 

We shall consider the specific limit ~> oo for investigating behavior of the doped fermions in the 
ordered phase. In this limit, the undoped groundstate \GS;c) = |G5(a, 6, c, c?))(c2^oo) satisfies 

Vc{x)\GS; c) = \GS- c), for aU x (48) 



and it is shown in Fig 6(a) Let us dope fermions in the 8v background. Fermion part of the 
Hamiltonian is the same with that of the IGT Eq.Q and the total Hamiltonian of the doped q8v 
model is given by 

n^nqsv + n^- (49) 

The constraint on the state H4U|) changes to 

al{x)al{x - i)al{x)al{x - 2) • e'''^^' ^'=^^^'=^\phys) = \phys). (50) 

In the ordered phase, the gauge variables cr](a;)'s take a definite value at each link. This means 
that the variables crf(x) fiuctuate very strongly, i.e., the ordered phase of the qSv model is nothing 
but the confinement phase of the IGT. This fact is pictorially observed for the limit oo. In 



this case, the undoped system is full of c-type vertices (see Fig 6(a) ). Then let us dope two fermions 
into the state \GS;c). In the lowest-energy state (i.e., the highest-weight state), the two fermions 
reside on the same site if their internal quantum number a are different whereas the fermions of the 
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same quantum number reside the nearest-neighbour (NN) site. In the former case, the gauge-field 
configuration is the same with the undoped state \GS; c). On the other hand in the latter case, new 
types of vertex appear at the sites occupied with the fermions as dictated by the local constraint 
(I50|l . Then energy of the resultant gauge configuration is higher than that of \GS; c) by an amount 
c"*, which is directly calculated from Hy. When one of fermions in a pair hops from its original 
position to a site with a distance L (on the lattice), L vertices change from the c-type vertices to 
the other ones because the hopping term of the fermion is accompanying the gauge operator <jf{x) 
which flips the arrow on the link {x,t). Then a separated pair of fermions cost energy ~ Lc^. This 
means nothing but the confinement of the doped fermions. See Fig^l 



(a) — > oo configura- 
tion of qSv model. 



(b) Doped fermions on 
the same site in the 
c'^ — > oo configuration. 



(c) Doped fermions in 
the ^ oo configura- 
tion. 



Figure 6: — > cxd configurations. 

Effective Hamiltonian of fcrmionic degrees of freedom in the ordered phase can be derived as in 
the previous discussion on the strongly-coupled IGT. For states with pairs of fermions residing the 
same sites like Ea. (|^ . 

Hf^{E) ^ Vc-HVc + VciiQcQciE - QcHQcy'QcnVc. (51) 

where (HqSv + 'H-4,)\tp) = E\iIj), Vc — Hi 'Pcix) and Qc — I ~ Vc- In the leading order of cP, 

= — ^[Pa/3(a;)P„/3(a; + l) + Pa/3(x)F„/3(a; + l)]. (52) 

^From (I52II and the discussion in Sec. 2, it is obvious that the superconducting phase of the fermions 
is realized in the ordered phase of the doped q8v model. On the other hand in the disordered phase. 
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no SSBs occur. We can expect that the phase boundary of the order-disorder phase transition of the 
background q8v system and that of the superconductivity of the doped fermions coincide, though 
the region of the deconfinement phase (i.e., the disordered phase with the translational symmetry) 
is enlarged by the fermion doping as in the EIGT. The phenomenon that the confinement induces 
SSB of internal symmetry is generaUy observed in various gauge theories. 

5.4 The limit = d'^ ^ oo 

It is interesting and also instructive to see how the doped fermions behave in the topologically 
ordered phase ///. To this end we consider the hmit = d? ^ co with keeping a^h — 1. We set 
c — d — K and take the hmit X — > 00. In this limit, the potential term Hy tends to 

Hv = K^Y.^Pa{x)+n{x)){Va{x + ai)+n{x + ai)) 

x{Va{x + a2)+Vb{x + a2)){ra{x + ai+a2)+rb{x + ai+a2)). (53) 

^From (|53() . the states of the lowest (vanishing) eigenvalue of i/y, \GS\cd), are specified by the 
projection operator Vcdix), 

Vcd{x)=Vc{x) + Pd{x), (54) 

Vcd{x)\GS; cd) = \GS\ cd), for all x. (55) 

There are various states, besides \GS;c) and \GS;d), which satisfy the condition H55|l . As the flip 
operator HfUp in (|41l) is applied on the \GS]cd)^ the a and/or b vertices appear. Therefore in the 
limit K — > 00, all various states in \GS;cd) are independent and degenerate. Among them, some 



symmetric configurations are depicted in Fig 7(a) and |7(b)^ i.e., the c-d columnar and c-d staggered 
configurations, respectively. 

Let us dope two fermions at site x in one of the various states and move one of the fermions from 
X to site y. The background 8v states change from the original configuration by the fermion hopping. 
For the various state in \GS;cd), one can verify that movement of doped fermions is restricted to 
straight lines, i.e., the horizontal and vertical lines. Otherwise a and/or b vertices appear as a result 
of fermion hopping (see FiglSJ. In this sense, fermions are deconfined only in one direction. In the 
present limit, the fermion system essentially reduces to the one-dimensional one, and therefore no 
SSBs take place as dictated by the Coleman theorem. This limit is in sharp contrast with Kitaev's 
point a = b = c = d= l. At Kitaev's point, there are no potential terms of the background 8v 
configuration and we can simply put tT|(a;) = 1 at all links. Then the doped fermions move without 
any interactions. 
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(b) Staggered 



Figure 7: = — *■ oo configurations of q8v model. 

6 Conclusion 

In this paper we have studied behavior of the doped fermions in the dimer and qSv models. In most 
of discussion we employed the IGT description of the models. We showed that the translational sym- 
metry breaking in the background dimer or 8v configuration induces SSBs of the internal symmetry 
in the doped-fermion sector. Similar phenomenon is well-known in the quantum chromodynamics 
in which confinement of quarks induces the SSB of the chiral symmetry of quarks. Result of the 
numerical studies on phase structure of the models, which shows that change of the background 
configuration really induces SSBs of the doped fermions, will be reported in the near future [T7|. 
It is also very interesting to construct more realistic models similar to the models in this paper 
for, e.g., strongly-correlated heavy fermions and compare phase structure of the models, etc. with 
experiments. 
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(c) Doped fermions in 
(c — d) columnar state 
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(b) Doped fermions in 
\GS;d) 



X 


4^ 














— »— z 


s— «— z 


i > L 


s— « — 




H 


H 





(d) Doped fermions in 
(c — d) staggered state 



Figure 8: Fermion hopping and change of background 8v configurations. Fermion hopping in one- 
dimensional straighthnes changes the c-vertices to the d-vertices. In the limit = ^ oo, resultant 
configurations of the 8v model have the same energy with the original ones. Therefore doped fermions 
move freely in the horizontal and vertical lines. 
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